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Abstract 

We consider the role of the perturbative Pomeron-to-parton splitting in 
the formation of the diffractive parton distributions. 

Diffractive deep-inelastic scattering (DDIS), 7*p — > X + p, is characterised by a large 
rapidity gap (LRG) between the cluster X of outgoing hadrons and the slightly deflected proton, 
understood to be due to 'Pomeron' exchange. Let the momenta of the incoming proton, the 
outgoing proton, and the photon be labelled p, p' and q respectively; see Fig. (Ha). Then the 
basic kinematic variables in DDIS are the photon virtuality, Q 2 = —q 2 , the Bjorken-x variable, 
xb = Q 2 /(2p • q), the squared momentum transfer, t = (p — p') 2 , the fraction of the proton's 
light-cone momentum transferred through the LRG, xp = 1 — p' + /p + , and f3 = xs/xp. 

It is common to perform analyses of DDIS data based on two levels of factorisation. First 
the i-integrated diffractive structure function may be written as the convolution of the 

usual coefficient functions C2 a as in DIS with diffractive parton distribution functions (DPDFs) 
a D [1]: 

F 2 D(3) (z P ,/3,Q 2 ) = J2 P f a B (w,z,t4), CD 

a=g,g J P Z \ z / 

with the factorisation scale fip usually taken to be Q. The DPDFs a D = zq D or zg D satisfy 
DGLAP evolution in \ip. The convolution variable z € [/3, 1] is the fraction of the Pomeron's 
light-cone momentum carried by the parton entering the hard subprocess. At leading-order (LO) 
the coefficient functions are C2, q (x) = e 2 5(1 — x) and C2 : g{x) = 0. The collinear factorisation 
theorem (Q]) applies when fipis made very large; it is correct up to power-suppressed corrections. 
In the second stage, Regge factorisation is usually assumed, such that the diffractive parton den- 
sities aP are written as the product of the Pomeron flux factor /p = Jdt exp(i?p t) x^~ 2ar ^ 
and the Pomeron parton densities a ¥ = zq ¥ or zg ¥ , that is, 

a D (xp,z,/i 2 7) = /p(xp) a ¥ (z,n 2 F ). (2) 

The Pomeron trajectory is ap(i) = ap(0) + apt. For simplicity, we omit the contribution of 
secondary Reggeons to the right-hand side of (f2]). Such an approach says nothing about the 
mechanism for diffraction: information about the diffractive exchange ('Pomeron') needs to be 
parameterised at the input scale ^0 an d fitted to the data. 
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Fig. 1 : (a) The perturbative resolved Pomeron contribution, which is the basis of the perturbative QCD approach, (b) 
the perturbative direct Pomeron contribution, and (c) the non-perturbative resolved Pomeron diagram, which accounts 
for the contribution from low scales, /i < /io- 

An alternative way to describe DDIS is to consider the heavy photon transition to qq or 
effective (qq)g dipoles, which then interact with the target proton via two-gluon exchange [2,3]; 
the case of the qq dipole is shown in Fig. [Hb). Here the Pomeron is modelled by the t-channel 
colour-singlet pair of gluons and we have an explicit form of the Pomeron parton distributions 
aF(z,fj?) given at the initial virtuality fj, 2 = k 2 /(l — z) fixed by the transverse momentum kt 
of the outgoing components of the dipole. Such an approach relies on the existence of a large 
saturation scale Q| > 1 GeV 2 [4] to act as an infrared cutoff and suppress the contribution from 

large dipole sizes, thereby justifying the use of perturbative QCD for the whole F^ 3 \ However, 
more sophisticated dipole models [5] now give a lower Q 2 S < 0.5 GeV 2 in the HERA kinematic 
regime, therefore significant non-perturbative contributions to inclusive DDIS are to be expected. 

The two theoretical frameworks are essentially contradictory: the Regge factorisation ap- 
proach is motivated by 'soft' Pomeron exchange, whereas the two-gluon exchange approach is 
motivated by perturbative QCD. Nevertheless, a 'doublethink' mentality exists whereby the two 
approaches are both commonly applied separately in the description of data, often within the 
same paper, with few attempts made to reconcile them. 

The two approaches can be combined [6-8] by generalising the 7* — > qq and 7* — > qqg 
transitions to an arbitrary number of parton emissions in the final state, as shown in the upper 
half of Fig. Q2 a). The perturbative Pomeron is described by a parton ladder ending in a pair of t- 
channel gluons or sea-quarks; the former is shown in the lower half of Fig. [Ha). The virtualities 
of the i-channel partons are strongly ordered as required by DGLAP evolution: /i 2 , <C . • • <C 
fj 2 <C . . . <C fip. The scale [i 2 at which the Pomeron-to-parton splitting occurs can vary between 
[1q ~ 1 GeV 2 and n 2 F . For fi < /xq> tne representation of the Pomeron as a perturbative parton 
ladder is no longer valid, and instead, in the lack of a precise theory of non-perturbative QCD, 
we appeal to Regge theory where the 'soft' Pomeron is treated as an effective Regge pole with 



intercept op(0) ~ 1.08; see Fig.QJc). 

The probability to rind an appropriate pair of t-channel gluons with transverse momentum 
It, integrated over l t , is given by the usual gluon distribution of the proton obtained from the 
global partem analyses, xpg(xp, fj, 2 ), where /U 2 = k 2 /(l — z) is the virtuality of the first t- 
channel partem in the upper part of the diagram. The emitted parton at the edge of the LRG in 
Fig- Etta) has rapidity r/ max and transverse momentum k t and carries a fraction (1 — z) of the 
Pomeron's light-cone momentum. For inclusive DDIS we must integrate over k t , accounting for 
the components of the Pomeron wave function of different sizes ~ 1/kt, which translates to an 
integral over ji of the form 



a D (x P ,z,Li 2 F ) = I -J- 



r> «s(M ) / ; 
R g x P g(x P ,iJ, 



fJ. Xp 



a p (z,^;/i 2 ). (3) 



The term /p(xp; /i 2 ) = [. . .] 2 /xp plays the role of the Pomeron flux in (0. R g is the skewed 
factor which accounts for the fact that in the lower parts of Fig. (Ha,b) we deal not with the 
diagonal but with the skewed (or generalised) parton distributions. For low values of xp 1 
this skewed factor is given by the Shuvaev prescription [9]. The notation a ¥ (z, [i 2 F ; fi 2 ) for the 
Pomeron parton densities means that they are DGLAP-evolved from an initial scale /i 2 up to the 
factorisation scale \j? F . 

At first sight, the integral ® appears to be concentrated in the infrared region of low \i. 
However, for DDIS we consider very small xp values. In this domain, the gluon distribution 
of the proton has a large anomalous dimension. Asymptotically, as xp — > 0, BFKL predicts 
xpg(xp, fi 2 ) ~ (/x 2 ) - 5 for fixed as [10]. In this case the integral ([3]> takes the logarithmic form, 
and we cannot neglect the contribution from large scales \i. 

By differentiating (O with respect to ln(/x|i), the evolution equation now reads [7] 



P a ,a>®a ,u + P a ,p(z)fp(xp;n z ). (4) 

a'=q,g P=G,S,GS 



Here the first term, involving the usual parton-to-parton splitting functions, P a>a ', arises from 
DGLAP evolution in the upper parts of Fig.[jla,c). The second (inhomogeneous) term, involving 
the Pomeron-to-parton splitting functions, P a ,p(z) = a p (z, /j 2 ; /j 2 ), arises from the transition 
from the two t-channel partons (that is, the Pomeron) to a single t-channel parton in Fig. QIa). 
The notation F = G,S, GS denotes whether the uppermost two t-channel partons are gluons 
(P = G) or sea-quarks (P = sfj] , while the interference term is denoted by P = GS. The 
LO Pomeron-to-parton splitting functions, P a ,p, were calculated in [7], where the perturbative 
Pomeron flux factors, /p(xp; fi 2 ), are also given. 

Simultaneously, we need to add in (Q~|) the direct Pomeron-photon interaction described 
by the coefficient function C2,p = Ct,p + Cl,f corresponding to the hard subprocess shown in 



1 Note that besides the two-gluon exchange we also account for the qq t-channel colour-singlet state, since even at 
very small xp the sea-quark contribution does not die out. Moreover, the skewed factor R q for the qq t-channel state 
is more than three times larger than that for the gluons [9]. 
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Fig. 2: The breakdown of the resolved Pomeron contribution for /i > 1 GeV to the total F 2 D ' 3 ' as a function of /x 2 
for two representative values of (xp, /3, Q 2 ). The integral over m(/i 2 ) is shown by the numbers in parentheses in the 
legend. Also shown are the integrated contributions from the the resolved Pomeron with fi < 1 GeV and the direct 
Pomeron with /i > 1 GeV. The secondary Reggeon contributions are negligible for the values of xp chosen here. 



Fig. [lib): 

if (3) = £ C 2 , a ® a D + Yl C 2,p- (5) 

a=q,g ¥=G,S,GS 

For the LO Ct,f the light-quark contributions in the limit /x 2 <C Q 2 are subtracted since they 
are already included in the first term of ([5]> via the inhomogeneous evolution of the DPDFs. This 
subtraction defines a choice of factorisation scheme. There is no such subtraction for the LO 
Cx p, which are purely higher-twist, or for the heavy quark contributions since we work in the 
fixed flavour-number scheme (FFNS). 

Thus, we see from © and Q that the diffractive structure function is analogous to the 
photon structure function, where there are both resolved and direct components and where the 
photon PDFs also satisfy an inhomogeneous evolution equation. 

As usual the evolution starts from some not-too-small scale /zo and all the contributions, 
both perturbative and non-perturbative, coming from fi < fio are parameterised in terms of the 
Regge factorisation as some input which should be fitted to the data. Note that inclusion of the 
inhomogeneous term in ((U) and the direct Pomeron coefficient function in ((5]) does not add any 
new free parameters to the description of the DPDFs. The LO Pomeron-to-parton splitting func- 
tions are known [7] and at next-to-leading order (NLO) they can be calculated unambiguously. 
The numerical results [8] presented below were obtained by fitting the HI LRG data [11] with 
Mx > 2 GeV and Q 2 > 8.5 GeV 2 , adding statistical and systematic experimental errors in 
quadrature, and taking the input distributions at /Xq = 2 GeV 2 , which gave a x 2 /d.o.f. of 0.84. 

From Fig. [2] we see that integrating the perturbative contributions to F^ 3 ^ starting from 
/i = 1 GeV we collect up to a third of the whole diffractive structure function F^ 3 \ Of course, 
the numerous corrections (higher order as corrections, power corrections, etc.) are not negligible 
at such low scales as /x ~ 1 GeV Nevertheless, this fact indicates that an important part of the 
diffractive partem densities comes from the relatively small size components of the 'Pomeron'. 
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Fig. 3: The breakdown of the contributions to the slope 9f 2 D ' 3 ' jd In Q 2 as a function of j3. Contributions apart from 
the usual DGLAP terms start becoming important for /3 > 0.4. The peak at (3 = Q 2 / (Q 2 + 4m 2 ) ~ 0.51 is due to 
the threshold for charm production from the direct Pomeron-photon interaction. 



That is why fitting the DDIS data in terms of the Regge factorisation one needs the effective 
Pomeron intercept ap(0) ~ 1.12 larger than that (ap(0) ~ 1.08) measured in the purely soft 
hadron-hadron scattering. 

In Fig.[3]we show the breakdown of the different contributions to the slope dF^ 3 ^ /din Q 2 . 
As it is seen, the inhomogeneous term in the evolution equation is fairly large and starts becoming 
important for j3 > 0.4. Note that the direct Pomeron contribution at large [5, mostly the twist-four 
F^ 3 \ only gives an important contribution to the total Kj for j3 > 0.9, but the contribution 
to the Q 2 slope starts becoming important at moderate (3. Therefore, it is not possible to avoid 
the presence of this contribution by simply excluding data points with j3 > 0.8 from the fit, as is 
done in the HI 2006 analysis [11]. In the analysis of inclusive DDIS data, the diffractive gluon 
density is mainly determined by the derivative dF^ 3 ^ /dlnQ 2 . The presence of these additional 
positive contributions to the Q 2 slope at large (3 apart from the usual DGLAP contribution means 
that a smaller gluon density is required for z > 0.4 compared to the HI 2006 Fit A performed in 
the usual Regge factorisation framework [11]; see Fig. HI 

In our analysis we account for the twist-four F^ contribution coming from the direct 
Pomeron-photon fusion, that is, from the last term in (f5]) corresponding to Fig. (Hb). This con- 
tribution, which goes to a constant value as ft — > 1, was calculated in [2,3] and turns out to be 
numerically appreciable. For the coupling to massless quarks, this contribution to (f5]) takes the 
form [7] 

= /^V h lf A f 2/n2 M*p; f ¥ M- (6) 

J^l y/l - 4pn 2 /Q 2 
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Fig. 4: The MRW 2006 pQCD DPDFs compared to those from the HI 2006 Fit A and the HI 2007 Jets fit. All the 
DPDFs shown here are normalised to My < 1.6 GeV by multiplying by a factor 1.23 relative to My = m p , where 
My is the mass of the proton dissociative system. 



Differentiating with respect to ln(Q 2 ) gives a positive term from the integrand evaluated at /x 2 = 
Q 2 /(4/3), analogous to the inhomogeneous term in the evolution equation, in addition to the 
negative term proportional to — 1/Q 2 obtained from taking the derivative of the integrand. 

A recent analysis [12] includes the twist-four F®^ contribution from all k t in addition 

to the usual Regge factorisation formulae of (Q]) and (O. It is found that the twist-four F^ 3 ^ 
term gives a negative contribution to the Q 2 slope, such that a larger gluon distribution at high 
z is needed than in the usual Regge factorisation approach. However, the contribution from 
fi 2 = k 2 / (1 — (3) < /Xq is already included in the input Pomeron PDFs taken at a scale /Xq = 1.5 
GeV 2 , so some double-counting is involved. Moreover, the DGLAP evolution is not accounted 
for in the dipole cross section, and the inhomogeneous term in the evolution equation of the 
DPDFs is neglected. The future measurement of F® [13] will provide an important check of 



the calculations of the twist-four component. 

The smaller diffractive gluon density at high z found in our analysis, compared to the HI 
2006 Fit A, is preferred by the data on inclusive diffractive dijet production, j*p — > dijet X' +p, 
where at LO the dijet system originates from the outgoing qq pair in Fig. (Ha,c) and the rest of 
the hadronic system X' originates from the other outgoing partons in Fig. (Ha,c). However, the 
data still tend to be slightly overestimated [14]. This can also be seen in Fig. 0] where the 'MRW 
2006' gluon at high z is smaller than the 'HI 2006 Fit A' gluon, but larger than the 'HI 2007 
Jets' gluon obtained from a combined fit, within the Regge factorisation framework, to inclusive 
DDIS and inclusive diffractive dijet data [15]. The \ 2 for the 190 inclusive DDIS points increases 
from 158 (HI 2006 Fit A) to 169 (HI 2007 Jets fit) on inclusion of the dijet data. Therefore, the 
gluon density determined indirectly from the inclusive DDIS data, under the assumption of pure 
DGLAP evolution, is different from the gluon density preferred by the dijet data. The apparent 
tension between the inclusive DDIS and dijet data in the Regge factorisation approach is partly 
alleviated by the inclusion of the perturbative Pomeron terms. 

Note that in the 'MRW 2006' analysis [8] the parton distributions of the proton, represent- 
ing the lower parton ladders in Fig. SIa,b), were taken from a NLO global fit. However, since 
the Pomeron-to-parton splitting is only calculated at LO, it may be more appropriate to take the 
inclusive parton distributions from a LO fit, where the gluon density is much larger at small x 
and does not take a valence-like form at low scales. In this case, the inhomogeneous term in the 
evolution equation, and the other direct Pomeron contributions, would be enhanced, leading to 
an even smaller diffractive gluon density at high z. Note that even if LO DGLAP evolution is 
used in the lower parton ladders of Fig. [Ha,b), NLO DGLAP evolution may still be used for the 
evolution of the Pomeron PDFs in the upper parts of Fig. [TJa,c). 

A more direct way to observe the perturbative Pomeron contribution to DDIS is to study 
the transverse momenta of secondaries in the 'Pomeron fragmentation' region, at the edge of the 
LRG, as in Ref. [16]. In contrast with the Regge factorisation of Fig. (He), where the transverse 
momentum distribution of the partons inside a Reggeon is assumed to have an exponential form 
with a rather low mean value of intrinsic k t , in the perturbative case of Fig. [TJa) the k t -distribution 
of the lowest jet in Figfjja) obeys a power law, given by the integrand of ©. Therefore we expect 
a larger k t of the secondaries with the long power-like tail. 

So far, measurements of inclusive diffractive dijets at HERA have primarily been made in 
the kinematic region of small = xb/xp where the resolved Pomeron mechanism of Fig. QIa,c) 
gives the dominant contribution, and the contribution from exclusive diffractive dijets, Fig. [TJb), 
is negligible. However, a first measurement has been made of dijets in DDIS with a cut on 
f3 > 0.45 in order to enhance the contribution from exclusive diffractive dijets, 7*p — > dijet + p 
[17]. Within the HERA kinematic domain the sea-quark component of the Pomeron is quite 
important (see Fig. [2]). This statement can be checked by observing the diffractive high-£/jn dijet 
distributions in the exclusive — > dijet + p process. At LO we expect the ratio of cross 
sections for high-E^ dijets produced via the two-gluon (gg) and qq t-channel exchange to be [7] 



(da£ p /dE T ) 



¥=G _ 



81 I" a(l - a)Q 2 R g x F g(x F , // 2 ) 
4 [E%, + q(1 -a)Q 2 R q x F S(x w ,n 2 ) 



(da^ p /dE T ) 



(7) 



P=5 



where a is the photon's light-cone momentum fraction carried by the high-E^ 1 jet, and the scale 
fi 2 = a(l — a)Q 2 + Ej,. Therefore, measurements of exclusive diffractive dijets at large Et in 
photoproduction would probe the presence of the qq t-channel exchange. However, it remains to 
be seen whether a measurement of exclusive diffractive dijets is feasible without imposing a cut 
on p. 

In summary, we have shown how to obtain universal diffractive partem densities aP which 
can be used in the description of different diffractive processes. We emphasise that the pertur- 
bative QCD contribution originated by the inhomogeneous (last) term in (@]) is not small and, 
starting from a relatively low scale hq = 1 GeV, may generate up to a third of the diffractive par- 
ton densities. Simultaneously, an account of this inhomogeneous term, and the twist-four 
contribution, leads to a lower diffractive gluon density at high z in comparison to that obtained 
using the Regge factorisation hypothesis. The presence of the perturbative (large scale) contri- 
butions reveals itself in larger transverse momenta (with a power-like tail) of the secondaries 
observed in the 'Pomeron fragmentation' region. The colour singlet sea-quark pair exchange is 
an important component of the perturbative Pomeron, which plays the dominant role in exclusive 
diffractive dijet production with E T 3> Q 2 . Parameterisations for the DPDFs and the diffractive 
structure functions are made publically availableJl 
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